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SYSTEM AND METHOD FOR UTILIZING SHAPE ANALYSIS 
TO ASSESS FETAL ABNORMALITY 

Cross Reference To Related Applications 

[0001] This application claims the benefit under 35 U.S.C. § 1 19(e) of U.S. Provisional 
Application No. 60/490,540, filed July 29, 2003 and U.S. Provisional Application No. 
60/493,442, filed August 8, 2003, both of which are hereby incorporated by reference as if 
repeated herein in their entirety. 

Background Of The Invention 

[0002] Prenatal screening methods are routinely employed to assess the likelihood of fetal 
abnormalities, commonly referred to as birth defects. For example, Down syndrome or Trisomy 
21 is the most common cause of severe learning disability and accounts for approximately one 
half of all chromosomal anomalies in live born children. 

[0003] Current methods to screen prenatally for trisomy 21 involve maternal serum testing for 
biochemical markers and/or ultrasound evaluation of biophysical markers. Maternal serum 
screening involves the quantitative analysis of biochemical markers and risk assessment based on 
likelihood ratios derived from the population distributions of affected and unaffected 
pregnancies. Ultrasound evaluation, however, has historically involved visual observation of a 
fetal image and deciding empirically whether the image looks "normal" or "abnormal" (for 
example, whether the cerebellum appears as a banana sign for open spina bifida). This approach 
requires extensive experience in the "art" of ultrasound and the interpretation is necessarily 
subjective. 

[0004] Accordingly, there is a need in the art for a system and method that adequately evaluates 
the morphological changes observed with birth defects during prenatal screening. 

Summary Of The Invention 

|0005] Embodiments of the present invention provide for utilizing shape analysis to assess fetal 
abnormality. According to one embodiment, coordinates of points identifying a shape in a fetal 
image are received, coefficients of one or more mathematical functions that describe the , 
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identified shape are determined, and the determined coefficients are utilized as markers to assess 
fetal abnormality. 

Brief Description Of The Drawings 

[0006] FIG. 1 is a flow chart that depicts a process for utilizing shape analysis to assess fetal 
abnormality in accordance with an embodiment of the present invention. 

|0007] FIG. 2 is a flow chart that depicts a process for utilizing shape analysis to assess fetal 
abnormality in accordance with an embodiment of the present invention. 

[0008] FIG. 3 is a block diagram that depicts a user computing device in accordance with an 
embodiment of the present invention. 

[00091 FIG. 4 is a block diagram that depicts a network architecture in accordance with an 
embodiment of the present invention. 

[0010] FIG. 5 is a flow chart that depicts a process for utilizing shape analysis of a fetal head to 
determine risk of fetal abnormality in accordance with an embodiment of the present invention. 

[0011] FIG. 6 is a screen shot that depicts outlining of a fetal head in accordance with an 
embodiment of the present invention. 

[0012] FIG. 7 is a screen shot that depicts outlining of a fetal head in accordance with an 
embodiment of the present invention. 

[0013] FIG. 8 is a flow chart that depicts a process for utilizing shape analysis of a fetal brow to 
determine risk of fetal abnormality in accordance with an embodiment of the present invention. 

10014] FIG. 9 is a screen shot that depicts outlining of a fetal brow in accordance with an 
embodiment of the present invention. 

[0015] FIG. 10 is a screen shot that depicts outlining of a fetal brow in accordance with an 
embodiment of the present invention. 
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[0016| FIG. 1 1 is a screen shot that depicts outlining of a fetal brow in accordance with an 
embodiment of the present invention. 

[0017| FIG. 12 is a screen shot that depicts outlining of a fetal brow in accordance with an 
embodiment of the present invention. 

Detailed Description 

OVERVIEW 

[0018J FIG. 1 depicts a process for utilizing shape analysis to assess fetal abnormality in 
accordance with an embodiment of the present invention. Upon receiving coordinates of points 
identifying a shape in a fetal image (step 100), the coordinates are used to determine coefficients 
of a mathematical function or functions that describe the identified shape (step 1 10). These 
coefficients are used as markers to assess fetal abnormality (step 120). 

[0019| As shown in FIG. 2, once the coefficient markers are determined (step 200), they may be 
used by themselves or with other markers to assess fetal abnormality (step 210). A fetal 
abnormality may be assessed by comparing a patient's coefficient markers to reference data of 
coefficient markers by conducting a statistical analysis. The reference data may contain 
unaffected patients and/or affected patients. The statistical comparison could result in a risk of 
fetal abnormality, a likelihood ratio for a fetal abnormality or an index value that could be 
considered within range or outside of range for a fetal abnormality. 

[0020| The use of multidimensional coordinates allows for the evaluation of a shape as a whole. 
In one embodiment, a statistical shape analysis involves the tracing of an outline around the part 
of a fetal image to be analyzed. The points that make up this curve are then analyzed to derive a 
function that best fits the individualized points around the outline. The coefficients of this 
function may be considered random variables and can be determined for each evaluated image. 
The coefficients may then be analyzed using multivariate statistics to determine if they are 
outliers compared to the normal population. A subject shape that has coefficients outside the 
normal ranges observed in control shapes would indicate that the subject shape was significantly 
different than expected. 
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[0021) Described below are several embodiments within which the present invention may be 
implemented. 

ARCHITECTURE 

[0022] FIGS. 3 and 4 illustrate the components of a basic computer and network architecture in 
accordance with an embodiment of the present invention. FIG. 3 depicts user computing device 
300, which may be an ultrasound machine (3-D, 4-D or color), MRI or CAT scan machine, 
fetoscopy machine, workstation, personal computer, handheld personal digital assistant ("PDA"), 
or any other type of microprocessor-based device. User computing device 300 may include a 
processor 310, input device 320, output device 330, storage device 340, client software 350, and 
communication device 360. 

[0023] Input device 320 may include a keyboard, mouse, pen-operated touch screen or monitor, 
voice-recognition device, or any other device that accepts input. Output device 330 may include 
a monitor, printer, disk drive, speakers, or any other device that provides output. 

[0024] Storage device 340 may include volatile and nonvolatile data storage, including one or 
more electrical, magnetic or optical memories such as a RAM, cache, hard drive, CD-ROM 
drive, tape drive or removable storage disk. Communication device 360 may include a modem, 
network interface card, or any other device capable of transmitting and receiving signals over a 
network. The components of user computing device 300 may be connected via an electrical bus 
or wirelessly. 

[0025] Client software 350 may be stored in storage device 340 and executed by processor 310, 
and may include, for example, imaging and analysis software that embodies the functionality of 
the present invention. 

[0026] FIG. 4 illustrates a network architecture in accordance with an embodiment of the present 
invention. The network architecture allows the imaging and analysis functionality of the present 
invention to be implemented on more than one user computing device 300. For example, in one 
embodiment user computing device 300 may be an ultrasound machine that performs all of the 
imaging and analysis functionality of the present invention. In another embodiment, user 
computing device 300a may be an ultrasound machine that performs the imaging functionality of 
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the present invention, and then transfers image or coordinate data over network 410 to server 420 
or user computing device 300b or 300c for analysis of the data. The analyzed data could further 
be transferred to another user computing device 300 belonging to the patient or another medical 
services provider for testing with others markers. 

[0027] Network link 415 may include telephone lines, DSL, cable networks, Tl or T3 lines, 
wireless network connections, or any other arrangement that implements the transmission and 
reception of network signals. Network 410 may include any type of interconnected 
communication system, and may implement any communications protocol, which may secured 
by any security protocol. 

[0028] Server 420 includes a processor and memory for executing program instructions, as well 
as a network interface, and may include a collection of servers. In one particular embodiment, 
server 420 may include a combination of servers such as an application server and a database 
server. Database 440 may represent a relational or object database, and may be accessed via 
server 420. 

[0029| User computing device 300 and server 420 may implement any operating system, such as 
Windows or UNIX. Client software 350 and server software 430 may be written in any 
programming language, such as ABAP, C, C++, Java or Visual Basic. 

ANALYSIS OF FETAL HEAD SHAPE EMBODIMENT 
[0030[ FIG. 5 provides an example embodiment of the present invention in which the shape of a 
fetal head is analyzed to assess fetal abnormality. Fetal abnormalities identifiable through the 
use of the present invention may include, among others, Down syndrome, Spina Bifida, Trisomy 
18, Trisomy 13, frontal bossing, unbalanced translocation, other chromosomal abnormalities, 
heart abnormalities and abnormalities of any major body organ, structural abnormalities and 
craniofacial abnormalities. 

|003l| According to this embodiment, a transverse view of the fetal head is obtained by 
ultrasound (e.g., UCD 300a), saved as a bit-map image and transferred to another computer (e.g., 
UCD 300b). Then, in step 500, a user (e.g., user 400b) employs digitizing software (e.g., client 
software 350), such as TPSDIG, DigitX, CalExcel, DSDigit, Digical, Windig or MacMorph, to 
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create an outline of the fetal head to be analyzed. (Digitizing software provides coordinate data 
when a user clicks on a particular point in a bit-map image.) In creating the outline, the user may 
identify two landmarks on the bit-map image along the OFD (ocipito-frontal diameter) axis so 
that the image may be aligned against a consistent axis and to allow for a uniform assessment of 
points on each image. 

|0032] As shown in FIG. 6, prior to placing the coordinates on the outline, a series of 20 lines, 
each the length of the OFD axis, may be generated utilizing software employing general 
algorithms. The first line overlays the OFD axis; the remaining lines are centered at the mid- 
point of the first line and then are rotated at 9 degrees from the previous line in a 
counterclockwise manner. Once this "fan" is drawn, the user may place points, via the digitizing 
software, on the outline of the fetal head where it crosses the equidistant lines of the fan. FIG. 7 
shows the resulting 40 points representing the outline. Each small circle connected in the loop 
represents a point on the image clicked by the user to identify the outline, and the two cross-hair 
symbols next to the number "2"s represent the OFD positioning landmarks. The use of the "fan" 
allows for more consistency in creating the outlines on different images. 

|0033] In step 5 1 0, the user-identified outline points are converted into data coordinates by the 
digitizing software. TABLE 1 represents the landmarks and outline points of FIG. 7 as stored in 
an output file. The output file lists the coordinate data in a tabular format, with the x-coordinate 
listed first followed by the y-coordinate. As indicated in TABLE 1, the two landmark 
coordinates are listed first, followed by the 40 point coordinates. 
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TABLE 1 



LM=2 

401.00000 304.00000 
180.00000 347.00000 
CURVES=1 
POINTS-40 

401.00000 302.00000 
399.00000 322.00000 
393.00000 337.00000 
386.00000 353.00000 
375.00000 365.00000 
365.00000 376.00000 
356.00000 386.00000 
345.00000 395.00000 
334.00000 400.00000 
320.00000 402.00000 
307.00000 402.00000 
294.00000 401.00000 
283.00000 400.00000 
271.00000 397.00000 
260.00000 396.00000 
247.00000 391.00000 
236.00000 387.00000 
220.00000 381.00000 
209.00000 373.00000 
194.00000 364.00000 
182.00000 349.00000 
178.00000 330.00000 
183.00000 313.00000 
189.00000 297.00000 
198.00000 283.00000 
207.00000 270.00000 
221.00000 261.00000 
236.00000 254.00000 
248.00000 248.00000 
260.00000 243.00000 
272.00000 240.00000 
287.00000 236.00000 
301.00000 236.00000 
315.00000 235.00000 
331.00000 238.00000 
346.00000 244.00000 
359.00000 252.00000 
371.00000 264.00000 
382.00000 275.00000 
391.00000 288.00000 
IMAGE=Imagel . bmp 
ID=1 

SCALE=0. 140234 
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(0034) In order to ensure the coordinate data is aligned properly with the two specified 
landmarks, the user may employ a software program that utilizes general algorithms to adjust the 
coordinate data by rotating and translating each data point so that the OFD axis lies along the 
horizontal axis and that the first point lies at the origin, (x,y) = (0,0). These adjustments do not 
change the shape of the outline being investigated. The software program may employ the 
following general algorithm: 

• Point 1 = Origin: Denote as Xi,Y| 

• Point 2 = Other end of OFD axis: Denote as X2,Y2 

• Determine Slope of Line: (Y 2 - Yi)/( X 2 - Xj) 

• Determine the angle of the slope: Theta = arctangent(Slope) 

• For each point in the shape rotate the point clockwise: 

• Subtract out the origin such that NewX = X - Xj, NewY = Y - Y\ 

• FinalX = NewX*cos(Theta) + NewY*sin(Theta) 

• FinalY = newY*cos(theta) - NewX*sin(theta) 

[0035] The set of FinalX, FinalY values represent the rotated data points such that the OFD axis 
lies horizontally. In different scenarios, it may be appropriate to rotate the image in a counter 
clockwise manner; if this is the case then the formulas for FinalX and FinalY are: 

• FinalX = NewX*cos(Theta) - NewY*sin(Theta) 

• FinalY = NewX* sin(Theta) + NewY*cos(Theta) 

[0036] In step 520, the user may employ software such as the NTSYSPC program (Exeter 
Software), which uses an elliptical Fourier analysis program to generate a series of best fit curves 
to the coordinate data using harmonics. The zero harmonic consists of only 2 coefficients and 
represents only translation (i.e. movement in the x and y direction) and not shape itself so these 2 
coefficients are usually not analyzed. Each harmonic comprises four coefficients of interpolation 
functions that describe the user-identified outline shape. In general, the more harmonics used, 
the better the fit to the coordinates; however, if too many harmonics are used the fit may be too 
good since small differences might appear due to poor placement of the outline point during step 
500. In this embodiment, the user evaluates the first harmonic (i.e., the four coefficients al, bl, 
cl and dl). 
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[0037] The user may choose alternatives to the method described above for aligning the 
coordinate data. For example, the user could rotate the image so that the long axis of the shape 
described by the first harmonic is parallel to the x-axis, or the user can align the image such that 
the first point of the outline is set equal to the end point of the first harmonic to the x-axis, or the 
user could combine both of these methods. These options are available as part of the Elliptical 
Fourier Analysis module in the NTSYSPC program. 

[0038] In step 530, the determined coefficients may be utilized as markers to assess fetal 
abnormality by conducting a statistical analysis to compare the patient's determined coefficients 
with reference parameters derived from reference data of coefficient markers (e.g., a statistical 
distribution of determined coefficients in the unaffected population and/or affected population). 
One exemplary method of doing this is to calculate the Mahalanobis Squared Distance (MSD) 
value for the patient's coefficients. 

[0039] Prior to determining the outline coefficients, the coordinates may be adjusted by scaling 
the coordinates so that the area enclosed inside the outline equals 1 . The original areas of the 
enclosed outline could be included as a separate variable in addition to the coefficients as part of 
the statistical comparison analysis. Also, many shapes change with growth in the fetus. As a 
result, it may be necessary to adjust the observed coefficients to account for gestational age of 
the fetus as part of the statistical analysis. 

[0040| Several other shape analysis methods besides elliptical Fourier analysis may be utilized 
by the present invention. All of the following methods determine coefficients of functions that 
may be used as markers to assess fetal abnormality. Some examples of outline methods are: 

A. polynomials 

B. cubic splines 

C. parametric polynomials 

D. parametric cubic splines 

E. bezier curves 

F. Fourier analysis of equally spaced radii 

G. dual axis Fourier analysis 

[0041] All of these methods deal with outlines, of which there are four classes - simple open 
outlines, simple closed outlines, complex open outlines, complex closed outlines. A simple open 
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outline is an open outline that has only one value of y for each x. A complex open outline is an 
outline that can have more than one y value for each x. Simple closed outlines are closed 
outlines such that if one draws a line from the center through the outline it only crosses the 
outline once whereas in a complex outline the line would cross the outline more than once. In 
most cases an open outline can be analyzed as if it is a closed outline by assuming there is a 
straight line from the last to first point of the outline or by mirroring all of the coordinates around 
the x (or y) axis, however the outline can be analyzed as an open outline. Some of the methods 
above, such as the polynomial method, work with open outlines while others like the elliptical 
Fourier analysis work with closed outlines. It is also possible to analyze 3-dimensional outlines 
(x,y,z) in accordance with an alternative embodiment of the present invention. 

[0042| To provide the reference data to which the determined coefficients are compared in step 
530, a statistical algorithm may be utilized to determine the statistical distribution of the 
coefficients in the unaffected population using a set of coefficients for a series of unaffected 
pregnancies. The distribution may be defined by a series of reference parameters for each 
coefficient or pair of coefficients such as means, standard deviations and correlations. The 
coefficients associated with an outline in a particular patient could be compared to these 
reference parameters and the chance that the coefficients could be equal to or more extreme than 
their observed values could be determined. 

[0043] If affected cases are available then they may be included in the reference data set to 
determine the statistical distribution of the coefficients in the affected population using a set of 
coefficients from a series of affected pregnancies. The patient's coefficients could then be 
compared to the distribution of coefficients in the affected population. Alternatively, the 
patient's coefficients could be compared to both the unaffected distribution and the affected 
distribution as defined by the reference parameters. For example, if both the unaffected cases 
and affected cases are multivariately normally distributed, the distribution function for the 
multivariate normal distribution can be used with the reference parameters for the unaffected 
distribution to determine a relative frequency for the unaffected distribution and then again for 
the affected distribution to determine a relative frequency for the affected distribution, and then a 
likelihood ratio may be determined. A likelihood ratio equals the quotient of the relative 
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frequency in the affected distribution to the relative frequency in the unaffected distribution. A 
risk result (e.g., 1 in 100) gives the chance that a patient with the same parameters could have a 
child with a fetal abnormality. A likelihood ratio gives the relative risk that the patient could 
have a child with a fetal abnormality. 

[0044] The likelihood ratio can be used to multiply a prior risk to determine a posterior risk 
(after accounting for the minor adjustment between odds and risk, if necessary). For example, 
the prior risk of Down syndrome is often based on maternal age. If the statistical distribution of 
the determined coefficients are independent of the distribution of other markers, the likelihood 
ratio could also be used to adjust the risk of Down syndrome determined by the other markers to 
determine the overall risk of Down syndrome based on the outline and the other markers. 
Alternatively, if the coefficients are not independent of the other markers, then risk of Down 
syndrome could be determined by utilizing reference parameters for a combination of the 
coefficients and the other markers together using multivariate normal distributions or other 
distribution functions. 

|0045| Examples of other markers include nuchal translucency, free Beta hCG and PAPP-A, 
Ductus Venosus, absent or hypoplastic nasal bone observed on ultrasound, maternal blood alpha- 
fetoprotein, maternal blood hCG, maternal blood unconjugated estriol, maternal blood dimeric 
inhibin A, maternal urine total estriol, maternal urine beta core fragment, maternal urine 
hyperglycosylated hCG, maternal blood hyperglycosylated hCG, ultrasound "soft markers" 
which include for example, nuchal edema or increased nuchal fold, short femur, hyperechogenic 
bowel, and echogenic foci in the heart, etc. 

[0046] As mentioned above, the other markers can be combined statistically with the results 
from the shape analysis to provide a final result to the patient. Alternatively, the medical tests 
for the other markers could be performed prior to the ultrasound exam and then, when the 
ultrasound exam is completed, the results of the other marker tests can be thereafter combined 
with the results from the shape analysis. 

|0047] For providing reference data according to this embodiment of the present invention, 
columns 2-5 of TABLE 2 A show the results of evaluating the first harmonic (i.e., the four 
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coefficients al, bl, cl and dl) of fetal head outlines in a study of 35 unaffected pregnancies in 
the first trimester. TABLE 2B lists the reference parameters consisting of the mean and standard 
deviation of each coefficient and the variance/covariance matrix consisting of the variance 
(standard deviation squared) and covariance between each pair of coefficients and the formula 
for calculating a Mahalanobis-Squared Distance (MSD) for each case. 
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TABLE 2A 
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TABLE 2B 



Variable | Mean Std. Dev. 
+ 

al | -.6189147 .0099677 

bl | .0104039 .0183979 

cl | -.0128969 .0130582 

dl | -.5127587 .00881 



Variance/Covariance Matrix (M) 





al 


bl 


cl 


dl 


al 


.000099 


-.000035 


.000064 


-.000088 


bl 


-.000035 


.000338 


-.000125 


.000036 


cl 


.000064 


-.000125 


.00071 


-.00006 


dl 


-.000088 


.000036 


-.00006 


.000078 



MSD = (X-jl) T M" 1 (X-|i) 

where : 

(X-p.) is a 4 element vector consisting of the patient's 4 
coefficients (al,bl, cl,dl minus their respective reference 
means) . 

(X~n) T is the transpose of the (X-[i) vector 
M" 1 is the Inverse of the 4x4 Variance/Covariance Matrix 
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[0048J TABLE 2C shows the MSD calculation in the 35 unaffected cases. A cut-off beyond 
95% of the observed data was established representing a point halfway between the last 2 MSD 
values (1 1 .714). TABLE 2D shows the results of 2 patients who happened to be carrying a fetus 
with Down syndrome, both of which based on their MSD calculation are outliers. As more data 
from Down syndrome pregnancies are gathered, additional reference parameters (e.g., means, 
standard deviations, and covariances) based on the Down syndrome cases could be calculated 
along with other statistical techniques such as likelihood ratios to determine the odds that a 
patient is carrying a fetus with Down syndrome. 
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TABLE 2C 



ID 




MSD 


1 


.5453869 


2 


.72770054 


3 


.86797516 


4 


.99110505 


5 


i 

X . 




6 


1 

X . 


1 7 S7 97 


7 


1 

X . 


"3 1 P 9 1 Q9 

J X O Z X j z 


o 


1 
X 


. D J 3 J / 4 


o 


1 

X . 


7 71 CM 4 Q 
/ / X J 4 4 J 


10 


1 

X . 


09 a A c 


\ 1 


o 

Z . 


JO / / O 




Z . 


u y U 41 0 / / 




9 

z . 


4 4 4 U DU 0 


14 


9 

Z . 


^"79^71 "7 
D 1 Z 0 / J. / 


1 J 


9 

z . 


D 4 Z / 4 Z / 




9 

Z . 


O Z / Oil" 


1 7 






1 8 

1 o 


r> 

j> . 


c f) 1 Q*7 £CL 
jUIO / DO 


10 




0 4 x 1 D / 0 


90 




7 fi P 77 Q6 
too 1 1 ? 0 


91 
i 


o 


y 1 0 zzuz 


99 


4 . 


0570078 


9T 


4 . 


5657224 


94 


4 . 


6067905 




4 . 


9834947 


zo 


5. 


5061446 


27 


5. 


5437792 


28 


5 


.877155 


29 


6. 


0931487 


30 


6. 


2047977 


31 


6. 


3209851 


32 


6. 


5548075 


33 


7 . 


4795723 


34 


9. 


0892238 


35 


14 


.338667 



TABLE 2D 



Outcome al_ bl. cj_ d! MSD 

36 .570043 -.015975 .01023 .558191 94.220565 

37 .58333 .076813 -.056937 .538497 403.86348 



- 16- 



PATENT 
Attny Dkt No: 51637-75 

|0049] In this embodiment of the present invention, only the data associated with TABLE 2B 
would have to be stored in the computing device for the statistical comparison analysis of a 
particular patient to be conducted. This would preserve storage resources in the event of 
reference data based on very large populations. 

ANALYSIS OF FETAL BROW EMBODIMENT 
10050] FIG. 8 provides an example embodiment of the present invention in which the shape of a 
fetal brow is analyzed using 3-D ultrasound to assess fetal abnormality. In traditional 2-D 
ultrasound the challenge to the sonographer is obtaining an image that is in the proper plane of 
view. The ability to consistently obtain the same angle and depth of view of the fetus requires 
subjective decision making during the ultrasound exam. Thus it is difficult to obtain consistent 
views of each fetus from exam to exam. This increases error when trying to analyze shapes from 
these images. 

|0051] 3D ultrasound allows for the simultaneous visualization of the fetus in 3 separate 2-D 
planes, as shown in FIG. 9. These three planes are called the sagittal (side-view), coronal (front- 
view) and transverse (top-view). In FIG. 9, the top left image represents the transverse plane, the 
top right represents the coronal plane, and the bottom right represents the sagittal plane. 

|0052] Misalignment of an image in any given plane distorts the images in the other two planes. 
Therefore, in order to perform proper 2-D image analysis of any given plane, the image should 
first be aligned properly in all three planes. For example, to analyze an image in the sagittal 
plane, a user can assure a proper view by aligning the fetus in the coronal and transverse views. 
In addition, by choosing a landmark in each plane (e.g. the white dot placed at the bridge of the 
nose in the FIG. 9) the depth of the section can also be defined. Using landmarks and aligning 
the fetus in all three planes can insure that the view of each fetus in any given plane is the same 
and thus reduce variation in shape due to differences in ultrasound technique. Ultimately this 
will improve the ability to see changes due to biological effects. 3D ultrasound further allows 
for image manipulation after the completion of the examination. Therefore, after capturing a 3-D 
image the operator can later rotate that image to the appropriate view. 
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[0053] The use of 3D sonography is a recent advance in prenatal ultrasound. The technique 
generates a multiplanar display of separate images in the coronal, sagittal and transverse planes 
obtained by the ascertainment of a single "volume". Once the volume is obtained the images in 
each plane may be rotated to provide consistent and reproducible planes as part of a "post 
processing" evaluation. This advance was previously unobtainable using conventional 2D 
sonographic techniques. Once a given desired image and plane is obtained, it may be 
superimposed on a digitized screen where geometric morphometric analysis may be performed. 

[0054] According to this embodiment, a fetal brow is analyzed in the sagittal plane - from the 
bridge of the nose to the midportion of the top of the skull. The fetal image of FIG. 9 is obtained 
by 3-D ultrasound (e.g., UCD 300a), saved as a bit-map image and transferred to another 
computer (e.g., UCD 300b). The sonographer (e.g., user 400a) places the landmark at the bridge 
of the nose in each plane. 

|0055] In step 800, a user (e.g., user 400b) employs digitizing software to create an outline of the 
fetal brow to be analyzed. As illustrated in FIG. 10, the user first places two landmark points on 
the fetal head in the sagittal plane - one at the bridge of the nose and the other at the midportion 
of the top of the skull. Once these landmarks are correctly in place, the user employs a software 
program that utilizes a general algorithm to find a point that has the same horizontal component 
as the landmark on the bridge of the nose, and the same vertical component as the landmark at 
the top of the head. Once this point is found, a series of 16 lines of equal length are drawn, as 
shown in FIG. 1 1 . 

[0056] The first line in FIG. 1 1 starts from the center point and goes through the landmark at the 
bridge of the nose. The remaining lines start at the center point and then are rotated at 6 degrees 
from the previous line in a counterclockwise manner until a line is drawn through the point at the 
top of the head which is at a 90 degree angle from the initial line. Once this "fan" is drawn, the 
user may place points on the outline of the skull where it crossed the equidistant lines of the 
"fan" using the digitizing software. This method ensures the user to get curves with less human 
error and more reproducibility. FIG. 12 shows the 16 resulting points placed in the sagittal 
plane. 
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[0057] In step 810, the user-identified outline points are converted into data coordinates by the 
digitizing software. The user may choose to use only the first eight points in subsequent analysis 
since these points more determinatively represent the fetal brow. 

[0058] In step 820, an elliptical Fourier analysis is employed using three harmonics (i.e., three 
sets of coefficients). The determined coefficients are utilized as markers to assess fetal 
abnormality by conducting a statistical analysis, such as a principal component (PC) analysis. 

[0059] In step 830, PC scores are determined based on a PC analysis of the determined Fourier 
coefficients, and in step 840, the MSD values of the PC scores are calculated. As discussed 
above, this shape analysis may also be combined with other markers to more completely assess 
fetal abnormality. 

[0060] The PC analysis utilized by this embodiment is a statistical technique used with 
multivariate data to reduce the number of variables used in further statistical analysis. (It can 
also be used as an exploratory analysis to see which of the variables are most important). The 
PC analysis is a standard statistical technique that generates a set of linear combinations of the 
underlying variables. These linear combinations represent new variables that can be used in 
other statistical analyses. The first linear combination is the most important variable and so on. 
Higher numbered principal components can be dropped from further analyses since they tend to 
represent noise. Once the principle components are determined (e.g., 0.3Xi + 0.2X 2 + 0.4X 3 ...), 
the PC scores can then be calculated which can be used as the variables in an MSD calculation. 

[0061] For providing reference data according to this embodiment of the present invention, 
TABLE 4 shows the results of evaluating the first three harmonics (i.e., three sets of the four 
coefficients A, B, C and D) of fetal brow outlines in a study of 32 normal pregnancies. 
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|0062] TABLE 5 shows the resulting PC analysis of the Fourier coefficients. 



TABLE 5 

PRINCIPAL COMPONENT ANALYSIS 

. pea al-d3 if outcome==0 & addnl_img ==0, mineigend .0) (obs=32) 
(principal components; 3 components retained) 
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Scoring Coefficients 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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-0 
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a3 t 


0 
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-0 
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-0 
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-0 
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0 
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0 
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0 
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-0. 
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0 
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-0 
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0. 


28823 


0 
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Variable I 
+ 


Mean 


Std. Dev. 


al I 


-1.53125 


.2485798 


bl I 


-.04715 


.0221123 


cl I 


.0037642 


.0151056 


dl I 


-.2098125 


.0330263 


a2 | 


-.000394 


.0018714 


b2 | 


.0091534 


.0040555 


c2 | 


.0765344 


.014202 


d2 | 


.0015418 


.0170104 


a3 j 


-.1625 


.0312441 


b3 | 


-.0149906 


.0066018 


c3 | 


.0026395 


.0092017 


d3 | 


-.0135719 


.0085003 
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|0063] In this embodiment, the first 3 principal components are retained and used to calculate 
the associated PC variables. First a z-score for each coefficient is calculated by subtracting the 
mean of the coefficient and then dividing the difference by the standard deviation of the 
coefficient from TABLE 5. Each PC variable is then created by multiplying the scoring 
coefficient for each (al-d3) by the z-score for each variable and then summing the 12 products, 
as shown in TABLE 6. 
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[0064] The PC variables are used in an MSD calculation to create an atypicality index, as shown 
in TABLE 7. The PC variables are independent of each other so the variance/covariance matrix 
is not needed. A cut-off of MSD> 7.81 ( 95 th percentile of the expected X 2 with 3 d.o.f.) was 
used to define outliers. 

TABLE 7 

Variable I Obs Mean Std- Dev- 

PCA1 I 32 -l.lSe-Ofl 2-441701 

PCA2 I 32 -2.71e-01 2.0337m 

PCA3 I 32 2.446-0=1 1. 01121b 

USD = (PCA1/2. 441701)a2+(PCA2/2.033714)a2+(PCA3/1. 01121b) a2 

|0065] TABLE 8A shows the coefficients from an elliptical Fourier analysis for 6 affected cases 
and 8 images that represent an additional image from a patient in the reference set. The MSD is 
calculated using the means and standard deviations of the coefficients, the PC scoring 
coefficients and the MSD formula for the 32 unaffected patients, the six affected patients and 8 
additional images from patients in the group. The results are shown in TABLE 8B. The 3 
outliers were all abnormal (trisomy 21, trisomy 18, and one fetus with multiple congenital 
anomalies). Two trisomy 21 and one case with an unbalanced translocation were not outliers. 
The affected cases were l-3x Trisomy 18, 2-21x2 campomelic dysplasia, 2-23x2, 2-23x3, 2-23x4 
(same patient) Trisomy 21, 865x2 Trisomy 21, 865-5x2 Trisomy 21, 865-17x2 translocation. 
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|0066| In this embodiment of the present invention, only the scoring coefficients, means and 
standard deviations of the coefficients from TABLE 5 and the standard deviations of the PC 
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variables in TABLE 7 would have to be stored in the computing device for the statistical 
comparison analysis of a particular patient to be conducted. This would preserve storage 
resources in the event of reference data based on very large populations. 

|0067| Several embodiments of the invention are specifically illustrated and/or described herein. 
However, it will be appreciated that modifications and variations of the invention are covered by 
the above teachings and within the purview of the appended claims without departing from the 
spirit and intended scope of the invention. 
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